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Abstract 

In this paper, we establish the full- wave electromagnetic scattering theory to study the 
electromagnetic scattering from infinitely long cylinders with radially anisotropic coatings. We 
show that the total effective scattering width can be dramatically reduced by the suitable 
adjustment of the dielectric anisotropy of the shell, while it is not the case for tuning the 
dielectric anisotropy of the core. Furthermore, we could make the cylindrical objects invisible 
when both dielectric and magnetic anisotropics are adjusted. In the long wavelength limit, 
we develop effective medium theory to derive the effective isotropic permittivity and perme- 
ability for the anisotropic coated cylinders, and the invisibility radius ratio derived from the 
full-wave theory for small coated cylinders can be well described within the effective medium theory. 
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I. INTRODUCTION 



In recent years, designing optical and electromagnetic (EM) invisibility cloaking devices 
has attracted great attention from both physics and engineering societies. Scientists have 
developed coordinate transformation method |l| and conformal mapping method ^ , which 
can protect the cloaked object of arbitrary shape from electromagnetic radiation. This idea 
has been verified by full-wave simulation of the cylindrical cloaking structure based on both 
finite element method (FEM) [3!| and finite-difference time-domain method (FDTD) {4]. 
The experimental demonstration at microwave frequencies ^ and even at optical frequen- 
cies ^1 for non-magnetic material has been realized lately. They have investigated the 
transformation based cloaks possessing spherical/cylindrical geometries with rotational sym- 
metry 

naaai 

|9| and geometries with reduced symmetries, such as eccentric elliptic 
cylinders lOj, square cloaks 11|], and arbitrary shapes 12|. In addition to the method of 
spatial transformation, Alu et. al. have proposed to use isotropic plasmonic coatings to 
render objects invisible [isl based on the dipolar cancellation, and proposed a parallel-plate 
metamaterial cloak [1^ that significantly reduces the total scattering cross section of a given 
two-dimensional dielectric obstacle in some frequency band. 

Due to the spatial compression, the required parameters should be anisotropic and de- 
pendent on the position [1]. To alleviate these contraints, EM invisibility cloaks have been 
realized by isotropic coatings based on the effective medium theory [l5|, [iGj. However, it 
still requires sufficient layers of alternating isotropic media in order to maintain the validity 
of that effective medium theory. Although recent development in metamaterial may allow 
the control of material parameters, such strict requirements for the parameters of a perfect 
cloak are still too difficult to realize in practice. As a way to reduce the number of required 
constitutive parameters to realize a manufacturable cloak, the incidence can be decomposed 



into transverse electric (TE^ 
with simplified parameters 



and transverse magnetic (TM) components so that 2D cloaks 
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18| can be used. Unfortunately, the simplified cylindrical 



cloak will cause large reflection at the cloak's outer boundary due to the impedance mis- 
match. From the discussion above, one can see that: (1) ideal cylindrical cloaks require 
6 parameters, some of which are infinite values; (2) simplified cylindrical cloaks requires 3 
parameters but still some value tends to the infinity at the inner boundary; (3) both ideal 
and simplified cylindrical cloaks need certain parameters to be inhomogeneous i.e., a func- 
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tion of position. This motivates our current work on achieving cyhndrical cloaks, which is 
based on only one single coating by a homogeneous and radially anisotropic medium. The 
invisibility condition can be established by the effective medium theory. It is obvious that 
the contraints in material fabrication are greatly alleviated. 

A lot of theoretical and numerical approaches have been developed to deal with the 
scattering problem of homogeneous anisotropic or even gyrotropic cylinders, e.g., integrodif- 



ferential equation 



volumetric integral equation method [2l|, combined field surface 



integral equation method 23|, finite difference method with measured equation of invariance 
(FD-MEI) 21[, and dyadic Green's functions 2^, |25|, |26|] . However, it should be noted that 



all these anisotropic tensors in those mentioned papers are defined in Cartesian coordinates, 
which is different from the radially anisotropy defined in the cylindrical coordinates as dis- 
cussed in this paper. Although the latter can be mapped into Cartesian anisotropy, then 
the anisotropic cylinder will be inhomogeneous and angle-dependent. 

Following our discussed motivation, in order to investigate the EM invisibility of radi- 
ally anisotropic coated cylinders, we develop the compact scattering theorem for a coated 
cylinder with homogeneous radial anisotropy, by extending the idea of embedding radial 
anisotropy in the orders of Bessel/Hankel functions in spherical case [27I. We focus our 
analysis on the invisibility characteristics of a homogeneous anisotropic coated cylinder. 
EM field components and scattering width of such coated cylinders are formulated. We 
discuss the roles of anisotropic parameters as well as the core-shell ratio on the reduction of 
the total scattering section. From the scattering algorithm, the effective permeability and 
permittivity of the core-shell system are also established so that the required condition for 
invisibility performance can be exactly determined. The numerical results are given. 



II. FULL- WAVE ELECTROMAGNETIC SCATTERING THEORY 

Let us consider the electromagnetic scattering from a radially-anisotropic coated cylinder 
of infinite length (see Fig. 1). 

For the simplicity, we assume a normally incident plane wave with the magnetic field 
along the z direction (i.e., the transverse magnetic (TM) wave). We assume that the coated 
cylinder is composed of the core with the radius a and the shell with the radius 6, and it is 
surrounded by free space (eo^yUo)- As for radial anisotropy, we indicate that the permeabil- 
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FIG. 1: Geometry of the scattering of a plane wave by a coated cylinder with permittivity and 
permeability of radially anisotropic tensors. The incident wave propagates along x-axis, and the 
magnetic field is along z-axis. 



ity and the permittivity tensors for the core and the shell can be expressed in cylindrical 



coordinates (f, 6, z) 



^ Spr ^ 



V 



Epe 
£ 



pz 



(la) 



( \lpr ^ ^ 



llpQ 
y ^ ^pz j 



(lb) 



where ep (or /x^) represents the permittivity (or permeability) tensor of the core by p = c 
and/or the shell by p = s. 



ikox—iu)t 



In Fig. 1, the incident magnetic field with unit amplitude can be written as Hz = e 
where ko = uj^Jeq^q. In what follows, e^*"^* will be omitted. In this connection, Maxwell 
equations in the core and shell are written as, 



V X H= - iuje ■ E 

V X EHu;^ ■ H. 



(2a) 

(2b) 
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In cylindrical coordinates, Eq. (2) can be decomposed into 



^^^prH^-^ - ^ (3a) 

lujfipeH^-^ - (3b) 

luiip^H^ — (3c) 

r or r o9 

-tuje ,eA ^^""^'^ - - ^ (3d) 
dr r 89 

-luSpeEff^— —. (3f) 



After some algebraic manipulations, we can obtain the governing equation of 

1 d OH 1 d^H, ^ rH -0 (4) 
Epg dr dr rSpr dO"^ ^ 

Inserting the solution in the form of Hz — ^(r)©(^) and introducing a seperation term m^, 
one obtains the ordinary differential equation for ^(r) 

For the part 0(6*), we have the form e'™^. The general solution to Eq. (5) is AJp^{kpr) + 
BNpm{kpr), where Jpm and Np^ respectively denote the mth-order Bessel and Neumann 
functions with the argument of kpr, kp — uj'^epo/jLpz, and {pmY — "^^ . Note that if it is 
an isotropic case (i.e., Spr — Spg), pm reduces to be an integer. 

In order to match the boundary conditions at cylindrical interfaces, the incident magnetic 
field can be expanded as follows, 

+ OC 

jjinc ^ ^ikox ^ ^ikor cose ^ ^ -m j^^j^^^ymO _ (g) 

m=— oo 

The scattered magnetic field in each region can thus be formulated 

+0O 

AmJcm{kcT)e''^^ T < a, (7a) 

m=— oo 

+ 00 

J2 I'^lBmJsmiksT) + C^iV,^(A;,r)]e^"^^ a < r < 6, (7b) 

m=—oo 

+ 00 

E ^"^[Mkor] + DmH';^\kor)y"'\ r > b, (7c) 



where Am, Bm, Cm, and Dm, are the unknown coefficients to be determined, and if, 



(1) 



represents the mth-order Hankel function of the ffist kind. 

Applying the boundary conditions of Eq and being continuous at r = a and r = b, we 
can derive the scattering coefficient Dm to compute the farfield parttern 

•Jcmikcd) JsmikgO) Ngm{ks(l) 

j^kJ^rrXkcd) -ksf^miha) -KN'^^ik^a) 

Jsm{ksb) Nsm{ksh) Jm{kob) 





Dr 



^ksJUksh) j^ksNUksh) koJmikob) 



(8) 



Jcm{kcCl) Jsm{kg(l) Ngm{kg(l) 

^kcJcmikcO-) ~ksJsm.{ksO) ~^sA^sm(^s'^) 

Js^iKh) Ngmihb) -H^\kob) 

jf-^ksJ'smiksb) j^kgNUKh) -k,Hl'\k,b) 

where the prime denotes the derivative with respect to the argument. Note that other 
coefficients can also be solved simultaneously, and hence the electric and magnetic fields in 
each region are obtained. The scattering problem for the transverse electric (TE) case, can 
be solved in a similar way and the corresponding scattering coefficients can be obtained by 
the duality of e — > /i and yU — > e. 



Scattering and extinction efficiencies are expressed through scattering amplitudes |28l |. 



oo 



o = A V ID I' 

^ sea , , / , -^m 



2 V ^ 

and Qext = tT / . ^e{Dm)- 
kob ^ 



m=—oo 



(9) 



m=—oo 



III. EFFECTIVE MEDIUM THEORY IN LONG WAVELENGTH LIMIT 



In this section, we present the formulation of our effective medium theory for the coated 
cylinders in the long wavelength limit, i.e., k^b <^ 1 and kgb <^ 1. As a result, the higher- 
order moments proportional to {kibY'-^^ {i = o, c, s) are expected to be neghgible and the 
effective scattering width of the coated cylinder is dominated by m = and m = 1 terms 
29l . ISOj. Thus, we set the conditions for an effective medium as Dq = and Di = 



29|, 
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corresponding to 



and 








-Jo{ksa) -No{ksa) 

-ksJoiksa) -ksN^iksa) 

Mk,b) No{k,b) Mkob) 

f^ksJ'^{ksh) f^ksK{k,h) koJ[,{kob) 



(10) 








0, 



(11) 



-JsiiksO) -Nsiik^a) 

-ksJ'si{ksa) -ksNi,i{ksa) 

Jsiiksb) Nsiiksb) Ji{kob) 

f-kj'^^{ksb) fLk^Ni,{ksb) koJ[{kob) 

where the subscripts si = \/e^eJe^ and cl = \fe^le^r ■ In the hmit of ki^b <^ 1 and k^b ^ 1, 
we can use the following approximations J<},{x) = 1, No{x) = ^hi(x/2), Jq{x) = — x/2, 
iV^(x) = I:, Ji{x) ^ f , J[{x) ^ i, Jlix) ^ ^, and Nl{x) ^ Substituting these 

approximations into Eqs. ( |TOl) and ( fTTi) and replacing eo, /^o by ^e//, /^e// respectively, one 
can obtain 



and 



Eseiicl ■ Ese + si ■ Ecd) 



(cl • Ese - Sl ■ Ece)] 



(12) 



(13) 



sl[(cl ■ E,e + sl ■ E,e) + (f )2-^i(cl ■ £se - sl ■ E,e)] ' 
It is evident that in the long wavelength limit, the electric and magnetic fields are decoupled. 
For an isotropic {p = c, s), they reduce to the results in Ref. [291]. As a 

consequence, the coated anisotropic cylinder in the long wavelength limit can be viewed as 
an effective homogeneous cylinder, and the scattering efficiency is expressed in a simpler 
form as follows, 



Qs 



TC 



An 



+ 2 



Eo - Eeff 



Eo + Eeff 



(14) 



Furthermore, the invisibility condition implies that Eeff and fieff in Eqs. f fT2|) and f|T3]) 
should be identical to the parameters of the host medium, i.e., Eq and /iq. Thus, the relation 
between core-shell ratio and the anisotropic parameters in the core and the shell can be 
drawn 



a 

r 



IJ's 



purely magnetic 



(15) 
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a 
~b 



(cl ■ Sse + si ■ ece){,ese - ■ si) 



l/(2sl) 



purely nonmagnetic. 



(16) 



(cl ■ Ese - si ■ ece){ese + ■ si) 
In fact, the invisibility conditions above possess a physical constraint, i.e., < ajh < 1 
according to the set-up in Fig. 1, which implies that given a certain set of anisotropic 
parameters for the shell and core, there may be no invisibility no matter how the core-shell 
ratio is tuned. It is also worth noting that if the coated cylinder is isotropic, Eqs. (fT5|) and 
f|T6|) respectively become 



a 

b 



/is - /ip 
/is - /ic ' 



(17) 



b [{es - ec)ies + eo) 
which are in accordance with the results in Ref. 3]. 



1/2 



(18) 



IV. NUMERICAL RESULTS 



Based on our theoretical results, we provide numerical calculations of the scattering 
efficiencies under different radial anisotropics and physical insights into the invisibility phe- 
nomena. 

Fig. 2 shows the full-wave scattering efficiency of a plasmonic cylinder coated with dielec- 
tric anisotropic shell versus the core-shell ratio a/6 at different sizes. First, for small coated 
cylinder, the long wavelength limit is valid, and hence one can resort to effective medium 
theory. In this connection, we discuss the effective permittivity (effective permeability is 
unity because of nonmagnetic components) in Fig. 2(a). We find that for certain radius 
ratios a/b, the effective permittivity equals Eq and the scattering cross section of the coated 
cylinder is almost zero correspondingly (see the dips in Fig. 2(b)). As a consequence, the 
coated cylinder is invisible or transparent. On the contrary, when the effective permittivity 
equals — eo in Fig- 2(a), the surface plasmon resonance takes place and strong scattering 
cross section is obtained (see the corresponding peaks in Fig. 2(b)). Then, we calculate 
the scattering efficiency for larger objects, which are shown in Fig. 2(c) and 2(d). In this 
situation, the long wavelength limit is not valid, and effective medium model cannot be 
used. In fact, with the object's size being increased, the scattering efficiency can still be 
reduced significantly (see the dotted lines in Fig. 2(c) and 2(d)) at the certain radius ratio 
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(we call "near- zero scattering" ratio), though these values are larger than those within the 
long wavelength limit. 




UH ■ 1 ' 1 ■ 1 ■ 1 ■ 1 -U.DH ■ 1 ■ 1 ■ 1 ' 1 ' 1 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 2: (a)The effective permittivity Eeff of the coated cylinder versus core-shell ratio for different 
Eso when b = O.OIAq. (b)-(d) The scattering efficiency versus core-shell ratio for different EsO and 
different values of b. Egg = 2.5eo (black dotted line); EgO = 5eo (red dashed line); Ege = 7.5eo (blue 
solid line), Other parameters are Ecr = £ce = — Seo, ^sr = 5eo, /^s^ = = /^o- 

Of particular interest is that significant reduction arises if we tune the value of Sge in the 
sheU. In Fig. 2(c), with decreasing (i.e., be reduced to 2.5£:o), the scattering efficiency 
can be reduced considerably, and thus the transparency or nearly "invisible" is attained for 
large objects. Fig. 2(d) have verified this tendency at an even larger size. In other words, 
through adjusting the dielectric anisotropy of the shell, it is helpful for us to realize much 
lower scattering cross section and better electromagnetic invisibility. Moreover, the near- 
zero scattering radius ratio can be tuned at the same time. We note that a resonant peak 
is also presented in Fig. 2, in the long wavelength limit (or for small size), the maximum 
of scattering is due to the resonance of D'^f^. In comparison with the small objects, the 
additional small peaks in Fig. 2(c) and (d) result from the contributions of D^2^ because 
the higher terms of the scattering coefficient cannot be neglected due to the increasing size 
of the objects. 
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FIG. 3: (a) The effective permittivity e^ff of the coated cyhnder versus core-shell ratio for different 
£ce when b = O.OIAq. (b)-(c) the scattering efficiency versus core-shell ratio for different values of b 
and different e^e- The values of See in (a)-(c) are: £cs = —Qeq (black dotted line); Ece = — 5eo (red 
dashed line); See = — 2eo (blue solid line), while other parameters in (a-c) are Scr = —5eo, e^e = 
Esr = 6^0 5 fJ-sz = fJ-cz = fJ-o- (d) the full- wave scattering efficiency versus core-shell ratio for 
b = 0.1 Ao with EcO = 2eo (black dotted line), EcO = 5eo (red dashed line), and EcO = 12eo (blue solid 
line), while the other parameters are Ecr = 5eo, EsO = £sr = — lOeo, ^J-sz = tJ-cz = /^o- 

Next, we consider a coated cylinder with radially anisotropic core but isotropic shell as 
shown in Fig. 3. Again, for a small coated cylinder, it is evident that the minimum of the 
scattering cross section takes place at the core-shell radius ratio determined by the condition 
Eeff = Bq. Hence the results based on the full-wave theory should be in accordance with 
those from effective medium theory. We can observed from Figs. 3(b-d) that the isotropic 
core results in the smallest scattering efficiency, i.e., better electromagnetic invisibility (e.g., 
red dashed curves in Figs. 3 representing isotropic cores in each case). This fact is further 
proved for larger sizes (see Fig. 3(c) and Fig. 3(d)), and even for either plasmonic core 
with dielectric shell (see Fig. 3(c)) or dielectric core with plasmonic shell (see Fig. 3(d)). It 
reveals that in the core-shell system incorprating radial anisotropy, the isotropy in the core 
is a better choice to minimize the scattering width, resulting in "good" invisibility. 
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FIG. 4: The full- wave scattering efficiency versus core-shell ratio for b = 0.2Ao. Parameters for the 
black dotted curve: See = £cr = — 3eoj e^e = 2.5eo> ^sr = 5eo, fJ-sz = fJ-cz = fJ-o- Parameters 
for the red solid curve Ece = Scr = — Seo, ^se = 2.5eo, ^sr = 5eO) Aisz = 0.5;Uo, /Xcz = — 7/io- 

It is known that for large coated objects, multipolar terms contribute to the scattering 
cross section, and the coated cylinder is visible for almost all the radius ratio (see Fig. (4)). In 
order to make the coated cylinder transparent or at least less visbile, one can further adjust 
the magnetic anisotropy to reduce the scattering cross section. From Fig. 4, we can conclude 
that the invisibility effectiveness of large objects can be greatly improved through suitable 
adjustment of both the dielectric and magnetic anisotropy and the scattering efficiency can 
be reduced for nearly one order in magnitude. 

If the size is sufficiently small, the core-shell ratio determined by the cloaking condition 
definitely results in better invisibility performance. In order to demonstrate the invisibility 
performance for a large size h = 0.2Ao improved by the adjustment of the magnetodielectric 
anisotropy, we present the patterns of magnetic field in the xy plane, calculated by finite- 
element solver of the Comsol Multiphysics. In Fig. 5(a), for the isotropic case (the same 
parameters as the red dashed line in Fig. 2(d)), there exists noticeable perturbation caused 
by the scattered field. Comparing Fig. 5(b) with 5(a), one can see that the perturbation is 
reduced through adjusting the anisotropic permittivity of the shell. If we take the perme- 
ability into account and properly tune its value, the near "invisibility" could be achieved for 
even relatively larger objects. In Fig. 5(c), the incident wave is almost unaltered outside the 
shelled cylinder, as if there was no scatterer. The specific core-shell ratios chosen in Fig. 5 
are based on the given material parameters, and obtained in numerical calculation for the 
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FIG. 5: Snapshots of the real part of total fields distribution in the xy plane around the coated 
cylinder with outer radius h = 0.2Ao at 0.3GHz. (a) Parameters are: Ece = £cr = — 3eo, £se = 
5eo, Esr = 5eO) fJ-sz = fJ-cz = fJ-o, a = 0.645&. (b)Parameters are: Eco = Ecr = — 3e0) ^se = 
2.560, ^sr = 5eoi Aisz = /^cz = fJ'O, o = 0.1556. (c) Parameters are: Ece = £cr = — 3eoi e^e = 
2.5eo, £sr = 5eo, tJ-sz = 0.5/io, fJ-cz = -7^o, a = 0.626. 



lowest scattering efficiency correspondingly instead of theoretical effective medium theory 
(due to the size constraint). It should be mentioned that the transparency mechanism in 
this paper is different from the ideal cloak proposed by Pendry l| via spatical compression, 
where there is no EM field in the core cylinder. Here the coated cylinder is penetrable. 



V. CONCLUSION 



In summary, we have demonstrated the full-wave scattering theory by magnetodielectric 
anisotropic coated cylinder of infinite length normally illuminated by a TM polarized plane 
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wave. The effective permittivity and permeability of the anisotropic coated cyhnder have 
also been derived in long wavelength limit. The non-scattering radius ratio obtained from 
the full-wave theory for small objects can be well described within the effective medium 
theory. Numerical results have shown that the effective scattering width c^ca = 26(5sca 
can be significantly reduced by adjusting the anisotropic permittivity of the shell, so as to 
achieve better transparency. However it is found that tuning the core anisotropy is not 
able to improve the efficiency of transparency. Furthermore, we could further improve the 
invisibility performance of the cylindrical objects by the adjustment of both electric and 
magnetic anisotropics. Although we have only considered the 2D case and losses materials, 
our research may be useful for the design of low-observability cylindrical targets. Based on 
the above investigation, we may extend our investigation to more general situations, such 
as obliquely incidence with arbitrary polarization, and finite-long cylinders. We could scale 
well towards the research of making multiple cylindrical objects "invisible" . 
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